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Given a suitably large and well connected (complex) graph state, any quantum algorithm can be implemented 
purely through local measurements on the individual qubits. Measurements can also be used to create the graph 
state: Path erasure techniques allow one to entangle multiple qubits by determining only global properties of 
the qubits. Here, this powerful approach is extended by demonstrating that even imperfect path erasure can 
produce the required graph states with high efficiency. By characterizing the degree of error in each path erasure 
attempt, one can subsume the resulting imperfect entanglement into an extended graph state formalism. The 
subsequent growth of the improper graph state can be guided, through a series of strategic decisions, in such a 
way as to bound the growth of the error and eventually yield a high-fidelity graph state. As an implementation 
of these techniques, we develop an analytic model for atom (or atom-like) qubits in mismatched cavities, under 
the double-heralding entanglement procedure of Barrett and Kok [Phys. Rev. A 71, 060310 (2005)]. Compared 
to straightforward postselection techniques our protocol offers a dramatic improvement in growing complex 
high-fidelity graph states. 



I. INTRODUCTION 

For certain algorithms, quantum computing offers the pos- 
sibility of exponential speed up over classical computing if 
significant obstacles to physical implementation can be over- 
come |T]. An important class of proposed implementations 
uses linear optical elements and photo-detection to perform 
the logical operations In these schemes the qubits are 
typically projected onto the required states using optical mea- 
surements. An undesired feature of this technique is that in 
general two-qubit gates are probabilistic. In a naive imple- 
mentation of the circuit model of quantum computation, this 
causes a decrease in success probability that scales exponen- 
tially in the number of two-qubit gates. It can be avoided by 
dividing the circuit into sub-routines, post-selecting success- 
ful implementations of sub-routines and then teleporting the 
sub-routine into the main algorithm llH 0,1111^. This approach 
is closely related to the one-way model of quantum comput- 
ing developed by Raussendorf and Briegel 0.[M@]> and has 
become a serious alternative to the circuit model. 

We divide the class of optical implementations of quan- 
tum computing into purely optical schemes where the logi- 
cal qubits are photons, and hybrid schemes where the log- 
ical qubits are matter systems. In the latter class, matter 
systems are used to store qubits, and optical excitations fol- 
lowed by projective photon detections implement the two 
qubit gates UJ, [U Qj, [H |l4, 15, 16, 17, 18J. Candidate 
matter qubits with optical transitions include NV centers in 
diamond f2^, quantum dots in microcavities fl^ or pho- 
tonic band-gap structures ||22ll . and neutral atoms in cavity 
QED I23l l24tl . An argument in favour of pursuing hybrid ap- 
proaches is that schemes using only linear optical elements 
require post-selection, and hence require a quantum memory 
||2|l- The natural candidates for such memories are typically 
matter systems, and should therefore implement the logical 
qubits. Optical delay lines are problematic as scalable quan- 
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turn memories, mainly due to absorption in the fibre, and low 
photo-detection efficiencies. Indeed, the problem of photon 
loss is a serious concern for any optical scheme, and the liter- 
ature contains a variety of proposals designed to circumvent 
this issue Ii25ll26ll . 

Although hybrid schemes do not suffer from qubit loss, they 
still have to solve the problem of photon loss during the en- 
tangling procedures. We can postselect successful outcomes, 
but photon loss can lead to a detector signature that, erro- 
neously, indicates a successful outcome. One hybrid scheme 
that achieves a degree of robustness against losing photons 
is the Double Heralding (DH) procedure proposed by Barrett 
and Kok 1 10]. The essence of this scheme is two consecutive 
parity measurements that side-step any photo-collection and 
detection inefficiency. In turn, this allows us to grow high- 
fidelity graph states for quantum computing. In this paper, we 
develop upon our previous work 1 19] to give a full model of a 
dominant source of errors in this scheme. We also extend our 
proposed adaptive growth strategy, which can significantly al- 
leviate the cost incurred by these errors. The solutions pro- 
posed here are suitable for correcting any parity measurement 
scheme where the error is known, and they move beyond mere 
post-selection where efficiency is traded for fidelity. 

We address the experimental challenge that optical schemes 
require indistinguishable (mode-matched) photon sources. 
The problem of mode matching can be divided into the match- 
ing of frequency, polarization, and spatio-temporal modes. 
These categories can be further subdivided into mismatching 
due to random fluctuations or due to fixed variations inher- 
ent in the configuration of the matter-qubit system. This latter 
category is especially important when the fabrication process 
for the nano-structure containing the matter qubit is not under 
complete experimental control (e.g., self-assembled quantum 
dots). It is well known that nano-structures generally possess 
a very broad distribution of properties that are an intrinsic con- 
sequence of the fabrication process, but that they can be accu- 
rately characterized or calibrated. Although here the details 
are worked out for temporal mismatch in the double heralding 
scheme, the strategies employed are more generally applica- 
ble. Indeed, any form of mode-mismatch can be tolerated, 
provided that the error is a known function of some detector 
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Bridging of large GHZ 
states into arbitary graph. 




FIG. 1: A schematic outline of the microcluster approach to graph 
state growth. Phase 1: GHZ states (or microclusters) are joined into 
larger GHZ states, in a manner that is inefficient for building GHZ 
states of arbitrary size. It yields a constant offline overhead to build- 
ing medium-sized GZH states. Phase 2: Once the GHZ states are suf- 
ficiently large, they can be merged efficiently into large GHZ states. 
Phase 3: The large GHZ resources can be used to grow any graph (in- 
cluding cubic lattices like cluster states, or the minimal graph state 
for a particular algorithm) via a bridging procedure. 



variable - that is, the error is monitored. 

More specifically, we consider the elimination of monitored 
errors that generate a so-called tilting in the ideally equally 
weighted amplitudes of a graph state. The error can be asso- 
ciated with a single graph vertex, and as such, can easily be 
represented in a generalized graph state notation that is intro- 
duced in section HTB] As illustrated in Fig.|2^, we assume that 
a large array of qubits is available. The device can perform 
parallel measurements between pairs of qubits, and the pair- 
ing of qubits can be optically switched. Following Nielsen |6tl, 
the proposed strategy proceeds by first constructing a resource 
of GHZ states, or microclusters: 
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and then fusing these into an arbitrary graph. This allows us 
to build a universal resource for quantum computing. 



Sections FlI Al fc HTBl will outline the problem with a brief de- 
scription of the double-heralding (DH) scheme, followed by a 
Jaynes-Cummings model [27] of cavity mismatch. Sections 
III Al fc HTBI discuss how to optimise the brute force growth of 
small, imperfect, GHZ states (see phase 2 in Fig.[T]i- The error 
under consideration will affect these states, such that through- 
out phase 1 , they have the more general form: 



\iji0a)) = COS{0)\O) 



sin(6')|l 



(2) 



Here 9 is the tilting angle of the vertex, and the vertex is un- 
tilted when 6 = n/A. Section III Dl shows that when con- 
structing a large graph with mismatched cavities, the ampli- 
tude distribution for A^-qubit entangled qubits does not dete- 
riorate past that of the distribution for 2-qubit construction — 
although there is a modest decrease in the intrinsic gate suc- 
cess probability. 

Section HlEl shows that when these GHZ resources are large 
enough to be used in phase 2 of graph state growth. Those that 
do not meet fidelity requirements can be purified probabilisti- 
cally. We call this removal of tilting errors, realignment. 

In section |lll] two procedures are described that can fuse 
the purified GHZ resources. With these procedures GHZ re- 
sources can be merged into larger GHZ states (see phase 2 in 
Fig. [Til, or bridged by a graph edge, allowing construction of a 
general graph (see phase 3 in Fig.[T]i- In the course of describ- 
ing these procedures it will be convenient to further generalize 
the graphical notation to include weighted graph edges and 
partial fusions. In section HVl we quantify the improvements 
gained over a naive postselection strategy. 



II. DOUBLE-HERALDING SCHEME WITH 
MISMATCHED PHOTON LEAKAGE RATES 

In this section we describe the double heralding scheme 
with both perfectly and imperfectly matched cavities. We 
construct a model for cavity photon leakage and discuss how 
this effects the double-heralding scheme. We discuss how the 
spread of the resulting tilting error can be limited, and eventu- 
ally removed by the realignment procedure. 



A. The double-heralding scheme 

The double-heralding scheme for the construction of graph 
states lfloll28ll uses matter qubits that have a three-level struc- 
ture denoted by the states |0), |1) and |e) [see Fig. ^p)]. The 
low-lying (|0), |1)) states are the computational basis states 
for the logical qubit, and |e) is an excited level that becomes 
occupied when a qubit in the |0) state is excited by a 7r-pulse. 
Of the two computational basis states, only the |0) is excited 
by the tt -pulse, as the transition from |1) is forbidden (e.g. by 
a selection rule ). 

The first stage of building a graph state requires an entan- 
gling operation between two qubits. The DH scheme for en- 
tangling two qubits consists of two rounds of measurements. 
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FIG. 2: (a) A schematic layout of a typical device that could implement our proposal. Illustrated here are an array of atom/cavity systems 
for which optical pumping and photon measurements can be performed in parallel. A form of optical switching will also be required so 
that measurements can be performed on arbitrary pairs of systems, (b) the energy level structure for the matter qubit, which allows only the 
transition |0) ^ |e) via 7r-pulses. (c) Photon leakage rates PA{t) and Psit) for two cavities A and B. If our scheme is not used, then to 
keep the errors within fault tolerance threshold of 1 — 10~^ the quantum computer must postselect on detector click times within a small time 
interval. 



and requires a resource of qubits in the |+) = (|0) + |l))/\/2 
state. If one could completely prevent photon loss, and have 
completely reliable number-resolving detectors, then only the 
first step of the DH scheme would be required — that is, 
the system would require only single heralding. For single 
heralding two matter qubits in different cavities A and B [see 
Fig. (|2^)] are prepared in the |+) state and then pumped with 
a TT-pulse. Upon relaxation the matter qubits will emit zero, 
one or two photons with probabilities 25%, 50% and 25%, re- 
spectively. Fig. illustrates how a beam splitter erases the 
path information of the photons. If one (and only one) photon 
is detected, and there is no photon loss, the system is pro- 
jected onto the |0, 1) ± 0) subspace, giving a successful 
entangling operation. 

When photon loss or lack of photon number information is 
included, a single detector click will project the system onto a 
mixture of the ideal result and the two photon subspace, |0, 0) . 
The unwanted part of the mixture can be eliminated by two 
more steps of the procedure. First, both matter qubits are ro- 
tated by the Pauli spin flip matrix X, which does not affect 
the desired part of the density matrix but converts the |0, 0) 
component into |1, 1). Next, the heralding is repeated by re- 
exciting the qubits with a 7r-pulse and waiting for any photon 
detection events. If this second round of heralding produces 
a single detector click then the procedure has succedded; the 
1 1, 1) component of the mixture is eliminated, and it is known 
that photon loss could not have occurred on the first round. 
This procedure relies on negligible dark counts in the photo- 
detectors ifToll . In order to avoid confusion, the term DH ap- 
plication denotes the two rounds of a successful entangling 
operation. 

After a successful application of DH the qubits are max- 
imally entangled in the state (|01) ± |10))/\/2- Applying 
the local rotation HX to either one of the qubits creates the 



graph state (|0+)±|1— ))/ \/2, where His the Hadamard gate. 
A larger graph can be constructed by an application of DH 
on two graph nodes that are both connected to only a single 
neighbour. If the DH application is successful, then the re- 
sulting graph is fused [28]. If the DH application fails, then 
the two qubits must be measured in the computational basis 
(a Pauli Z measurement) in order to remove them grom the 
graph state. More generally, if the qubits have more than one 
neighbour each, a successful DH application will yield a graph 
with one qubit having all the connections of the original two 
qubits. The other qubit (up to a Hadamard gate) will be con- 
nected to this qubit as a single dangling bond, or a "cherry". 



B. The effect of cavity mismatch 

The effect of mismatched cavity leakage rates on the 
double-heralding scheme will now be described in the Jaynes- 
Cummings model. We show that mismatched cavities cause a 
tilting in the amplitude of the resulting state, which depends 
on the photon detection times. With good time resolving 
photon detectors and calibrated photon sources this tilting is 
known to the experimenter and can be coiTected by the scheme 
presented here. In the non-ideal case, if atom A is in a cav- 
ity that tends to emit photons faster than the cavity that hosts 
atom B, then an earlier detector click means that the photon 
is more likely to have come from cavity A, hence achieving 
only partial path erasure. In order to know the resulting state, 
each atom-cavity system must be calibrated by measuring the 
photon leakage rate from a qubit prepared in the |0) state. We 
denote this probability distribution by Px{t), where the x in- 
dexes the atom/cavity system. Two examples of mis-matched 
PAit) and Psit) are shown in Fig. where perfect path 
erasure occurs only when the curves cross. 
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A suitably general model for the probability distributions 
Pj. (t) can be constructed by considering the three-level atomic 
system and its coupling to an electromagnetic cavity mode. 
With only one quanta of energy available to the cavity mode, 
its quantum state will be described in the Fock basis of |0) 
and 1 0), with no photons, and one photon of energy fkj, re- 
spectively. If we monitor a joint atom-cavity system x for the 
emission of photons, and none are detected, then its evolution 
can be described by a non-Hermitian conditional Hamiltonian 

^4 = g.(|e)(0|a+ |0)(e|at) _ Ijtj^, (3) 

where is known as the quantumjump operator, 

Jx — \J ■ (4) 

The constant represents the Jaynes-Cummings coupling 
strength between the cavity mode and the optical transition, 
and K,^ quantifies the leakage rate of the cavity (fi, = 1). Note 
that the last non-Hermitian term is responsible for the irre- 
versible evolution of the system. This irreversible decay also 
reduces the norm of the wavefunction N = (vE'|vE') at a rate 
TV that represents the probability of detecting a photon. In 
general, any measurement event associated with a jump oper- 
ator J will occur with probability JV|\1/). Defining the 
amplitudes of a system as follows: 

1*)^ = ci|e)^|0) +C2|O)^at|0) +C3|1)^|0), (5) 

the detection probability satisfies, N — k,x\c-i^ . We define 
Pxii-i Kx) as the solution for N{t) when Eq. ^ is solved 
with the initial condition ci{t — 0) = 1. Note that, since ci is 
coupled to the decaying component C2, it too will vanish over 
time. 

When a photon is detected the system is projected onto a 
new state by the quantumjump operator J^- In the case of 
a single photon source the constant factor has no effect, 
but when we model the double heralding scheme it becomes 
physically important. 

If we consider two of these systems A and B, then whilst 
no photons are detected the joint system will evolve according 
to J^^AB — J^A + '^B- If detectors Z?+ and £)_ are placed be- 
hind a 50/50 beam-splitter then the corresponding jump oper- 
ators J+ and J_ with each be a mixture of the photon modes a 
and b. To retain the same conditional Hamiltonian we require 
that; 

JaJa + jUb = jU+ + J-J- (6) 

This requirement alone does not uniquely define the new 
jump operators. With the beam splitter transformations in 
mind, it may seem nature to presume that J+ and J_ will 
be y^KX72(a + 6) and ^ nsj'^.ia — b). However, both detec- 
tors — being behind a 50/50 beam splitter — have an equal 
chance of registering a photon, hence we require: 

(*|J|J+|*) = (*|j1j_|*>, (7) 



which is not satisfied by the our first guess. Jointly, the two 
conditions in Eqs. (|6]l and (|7]i restrict the jump operators to: 

4 = (f)^±exp(^<^)(^)% (8) 

where the undetermined phase, exp(i0), is related to any dif- 
ference in path length. This phase is unimportant as it van- 
ishes after two rounds of double heralding. We therefore set it 
to 1 from now on. 

Using this formalism we can now turn to modelling the ef- 
fect of mismatch on the double heralding scheme. We first 
consider the preparation of qubits in the state |+)a|+)b and 
then apply the double heralding scheme. To pass the first 
round the detector must click at time ti, then be left without 
being interrupted by a second detector click, until the decay- 
ing amplitudes, Ai, A2, Bi, B2, are negligible. This is will 
generate the state: 

[PA{h)i\0)A\l)B+PBitl)^l)A\0)B) |0) 
^^^^"^ ^ VPA{h)+PB{h) ' 

(9) 

which lends itself to an intuitive interpretation. Each ampli- 
tude squared is simply the relative probability of one cavity 
having emitted a photon compared to either cavity emitting a 
photon. Proceeding with the the double heralding procedure, 
both qubits are flipped with XaXb, and another tt -pulse is 
applied. Upon detection of a second photon at time t2, the 
system is projected onto: 

\^fiep))AB =CO8{e„)\l)A\O)B+Sm{0p)\O)A\l)B, (10) 

where the amplitudes are represented as a function of an angle, 
6f^, which we call the tilting angle: 

s , PAitl)PB{t2) Y^ ,,,, 

When 9p — ±7r/4, the state simplifies to a standard graph 
state, otherwise we say the state belongs to a class of general- 
ized graphs that has a tilted vertex. As shown in Fig. ([3^), we 
represent the states as graph states but with the tilted vertex 
labelled by an angle 6. In a similar fashion to the construc- 
tive definition of pure graph states, we constructively define 
the quantum state associated with a tilted graph. Whereas a 
pure graph has every qubit initialised in the state |+), a tilted 
vertex is prepared in the state: 

16*) = cos(6l)|0) +sin(6l)|l) . (12) 

After this different preparation stage, control-Z gates are ap- 
plied to all qubits connected by a graph edge. 

Having demonstrated the modelling techniques used on a 
simple two qubit system, and having defined the graphical 
representation of a tilted vertex, we shall now state the results 
for more complex systems. 

Consider the effect of the double heralding scheme when 
the qubits used correspond to tilted vertices in tilted GHZ 




FIG. 3: The effect of a successful or failed application of the double heralding scheme (where the qubits used are within the dotted ellipse) 
for the non-ideal case: (a) when the initial state is two separable qubits in the untilted |+); (b) when the qubits are part of a tilted GHZ state; 
(c) when the qubits are part of a tilted GHZ state, but have had a Hadamard removed. All successful outcomes have a single tilted vertex 
parametrised by Oft, which is specified by Eqn. l ll6t . The probability of success is given by Eqn. l ll4t . (d) some identities for tilted vertices: 
when 6 = ±7r/4, the graph is equivalent to a pure graph; typically in this paper we only label Hadamard rotations as they have an important 
effect on the double heralding scheme, although in section HTdI we use the second identity illustrated — where an X rotation transforms a 
tilted angle 6 ^ -r /2 - 6. 



states. Fig. ^p) represents this in the new graphical nota- 
tion (where local rotations corresponding to detector parity 
have been ignored). From our constructive definition of tilted 
graph states it follows that the two initial states, a and 6, have 
the form: 

=cos(0x)|O)®"^ +sin(f?,)|l)f"^ (13) 

where x = a,b and Ua = ni, = n. Performing a parity mea- 
surement between a qubit in a with a qubit in b has a success 
probability: 

Pi9a, Ob) = coa^iOa) sin^Ob) + sin^Oa) cos2(0f,), (14) 

which is the upper bound of success for the double heralding 
scheme — reduced upon consideration of photon loss. Upon 
failure the two qubits are measured in the Z basis, projecting 
all 2n qubits into a separable state. Success yields a single 2n 
qubit entangled state: 

i^fm.y = cos(0^)io)rii)f" (15) 

+ sin((?^)|l)nO)f , 

which is, up to a spin flips on n qubits, equivalent to another 
tilted GHZ state with a tilted vertex parametrised hy 9p, such 
that: 

fn ^ . ts.n^{eb)PA{ti)PB{t2) Y'^ .... 



This can be shown to be consistent with the less general 
result when the initial states are pure graph states, by substi- 
tuting in da = Ob = TT /A and deriving (fTTl i. Furthermore this 
simplification is possible when da = 9b, which is a point that 
we shall return to in section (IIIDb . 

In Fig. Qc), the graphical notation is used to describe how 
the double heralding scheme applies when the qubits used are 
not the nodes of a GHZ state. Note that the notion of cen- 
tral node requires some clarification, since the state described 
by Eqn (fT3T l has no single qubit playing a privileged role that 
would make it 'the' node. Indeed, applying the double herald- 
ing scheme will have the same effect independently of which 
qubit is used. However, in Fig. (|3]c) the qubits used possess an 
additional Hadamard rotation compared to Eqn (fTsT l. This has 
an effect on the shape of the graph for successful and failed 
outcomes. However, it has no effect on the probability of suc- 
cess, or the value of the resulting tilted vertex, which both still 
obey Eqs. ( fT4l l and (fT6] l. 

So far we have given a full account of the effect of cavity 
mismatch when the detectors click at times ti and t2- What re- 
mains to complete the description is to specify the probability 
distribution for the random variables ti, t2- The probability 
density for a single detector click in round one is: 

Qi{h) = cos2(0„)sin2(0fc)P^(ti) (17) 
+ siii2(0,)cos2(0fc)PB(ii). 
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FIG. 4: Plots of E{F'^), the expectation value for the fidelity squared. The example cavities used are described in appendix. JB). (a) A 3D-plot 
of the solution for E{F^) as a function of sin^ [9a) and sin^ (&b). The red and green lines correspond to the cross sections used in (b). (b) Two 
cross sections of the 3D-plot, along sm^{9a) ~ sin^(^i,) (red) and sin^(6la) = 1 — sin^(Sj,) (green). 



Given an event at ti, the second round is governed by the 
probability density distribution: 



Q2{t2\tl) 



(X + Y)PA{t2)PB{t2) 
XPA{t2)+YPB{t2) ■ 



(18) 



where 



X = cos'' {0a) sin' {9b)PA(h)PB{t2) and (19) 



Y 



sm'i9a)cos'{9b)PB{ti)PA{t2), 



(20) 



These are convenient variables that will allow us to cast many 
expressions in a concise form. The product of these distribu- 
tions gives the joint probability density distribution for ti and 

t2- 

Qi2{.hM) = Qi{ti).Q2(t2\h) ^x + Y . (21) 



C. Expectation values for some measures of gate quality 

In this section we consider the expected quality of an at- 
tempt at double heralding. The results of this section will 
form the criteria for evaluating different strategies for graph 
growth that are discussed in the subsequent section. We con- 
sider two tilted vertices 6a and 6b, with two leakage rates 
PA{t) and PB{t), from which we can determine the proba- 
bility distribution for 6/3. Our quality measure must incor- 
porate both the probability of success, defined in Eqn. ( fT4l l. 
and some function that measures how close dp is to the ideal 
values of ±7r/4. Without loss of generality, we shall only con- 
sider +7r/4 as the ideal case, and assume known phase errors 
are always corrected with a local Z rotation. 



Although we could use a pure state entanglement measure, 
such as the von-Neuman entropy, we shall instead consider 
functions of the inner product between the tilted state and the 
ideal state: 



f{6p) = |(vl/(0^)|vl/(7r/4)>|^ 



(22) 



M^Sfj)). 



In this section we will give some plausibility arguments for 
the value of these measures. However they are ultimately es- 
tablished by the nature of the realignment procedure, which 
distils the tilted states and is the topic of section. HIE! 

Since even a failed double heralding application generates 
|5'(0)), which has /(O) = 1/2, we choose the quantity of 
interest to be F{6f3) = f{6fj) - 1/2. Using Eqn. (O, Op is 
eliminated from F{dp) to give: 



F : 



(23) 



X + Y' 

where X and Y have been defined in Eqn. ( fT9l ). To get the 
average fidelity over all possible detector click times for ti 
and t2, we multiply by (312(^1,^2), which takes the simple 
form Q12 = X + Y, and integrate over <i and t2'. 



E{F) = isin(20,)sin(206) 



[PA{t)PB{t)Y dt, 



(24) 

Hence the expected fidelity can be split into two independent 
factors, one dependent on initial tilting angles, and one depen- 
dent on the overlap between the probability distributions. 

However, many different fidelity distributions may have the 
same expected fidelity and yet represent very different re- 
sources. Consider two distributions: the set [/ of 2 n-qubit 
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tilted graphs all with sm{29u) = A/2; and the set V of 4 
n-qubit tilted vertices with sin(26'y) = A/A. If we sum the 
values of F they both equal A, but the two distributions are 
of different utility in subsequent attempts at double heralding. 
Attempting to double herald set U, will produce a distribution 
of 2n-qubit tilted graphs with expected fidelity proportional to 
sm^{29(j) — A? /A. Whereas the same calculation for set V , 
gives 2%\Y?(2Qv) = Heuristically, this indicates that 

quality is better than quantity. Quantitatively, it tells us that 
may be a better measure of gate quality. As alluded to 
earlier, this is a point that is supported by section. Ill El The 
expected value of F"^ is: 

^ ' 4 7 eiPA{tl)PB{t2) + Q2PB{tl)PA{t2) ' 

(25) 

where, 

Oi = cos2(0,)sin2(0f,), (26) 
92 = sin' {ea)cos^ {9b). 

Unfortunately, the general case does not admit tilting angles to 
be factored outside the integral. Furthermore, even for speci- 
fied cavities, evaluating the integral can be quite involved. We 
proceed by considering an example, for a particular pair of 
cavities. The details of these cavities are given in appendixiBl 
and are the same as used to generate the example probability 
distributions given in Fig. (|2]i. In appendix [A) we demonstrate 
the details of how to evaluate E{F''), and we give a general 
proof that for any probability distribution the behaviour of 
E{F'') is qualitatively the same. Hence, we can draw con- 
clusions from Fig. (|4|i without concern that they do not carry 
across for other pairs of cavities. 

D. Strategies for Phase one GHZ growth 

This section proposes strategies for the "phase 1" growth 
of GHZ states on the basis of as a measure of the ex- 

pected gate quality. Firstly, we consider the procedure given 
two tilted GHZ states with tilting angles {6a, 0^) and cavity 
leakage rates PA{t) and Pb {t). Secondly, we discuss the more 
complex issue of how best to divide 2N GHZ states from a set 
(6*1, ....02n) into N pairs of GHZ states — where the objec- 
tive is that the selection maximises E{F'') summed over all 
pairs. 

Consider the first problem, for which the only freedom we 
have that can change the outcome is the application of local 
rotations prior to a double heralding application. If we ap- 
ply an X rotation to one of the qubits, then up to some lo- 
cal rotations on other qubits in the graph this interchanges 
the magnitude of the qubits |0) and |1) — the tilting angle 
transforms 9 -k /2 — 9. Hence, if 9a and 9s are of simi- 
lar magnitude, they become dissimilar, and vice-verse. This 
also inverts the success and failure probabilities, so applying 
spin flips alternates between high and low success probabili- 
ties. When the success probability is high, one qubit is much 
more likely to emit a photon than the other, and we can expect 
the resulting graph to be more tilted. Interestingly, the X flips 



have no effect on E{F), as these two effects cancel out ex- 
actly. However, E{F'') is more sensitive to fidelity than suc- 
cess probability. Guided by this measure, we prescribe that a 
spin flip is applied when the tilting angle are initially far apart, 
I sin^(6'a) — sin^(6'h)| > 1/2. To put this argument on a more 
quantitative grounds, we consider when the tilting angles are 
symmetric about 7r/4, and hence far apart. Now is 
the green curve of Fig. (|4|i, and an X flip changes makes the 
tilting angles identical and changes E{F') to the red curve. 

We now turn to the second question: how best to pair up a 
set of tilted GHZ states (9a, 9),, ...92n) for double heralding. 
The optimal strategy is computationally hard because it re- 
quires evaluating a complex measure of success, like 
for every possible combination of pairs. Furthermore, since 
the quantum computer will have to perform these assessments 
whilst running, time is a critical factor, as delays between 
rounds of double heralding will increase the amount of de- 
coherence. Non-optimal, yet good strategies that are com- 
putationally efficient to implement may be divided into two 
categories: those based solely on matching cavities; and those 
based solely on tilting vertices. It may be the case that there 
also exist computationally efficient strategies that use both 
pieces of information, although whether this is the case is not 
clear at this point. 

In this paper we propose a strategy that takes only tilting 
angles into consideration. We have chosen this over a cavity 
based strategy for two reasons: (1) Although we may match 
cavities during GHZ growth, to make a unified graph we will 
eventually have to connect qubits from mis-matched cavities; 
(2) The problem does not admit a general answer, as it de- 
pends on the connectivity of the target graph, and the distribu- 
tion of cavities used. 

Having motivated our proposed strategy, its enunciation is 
simple. The tilted GHZ states are sorted into an ordered list 
from descending to ascending tilting angle, and the adjacent 
items on the list are paired up. Again this utilizes the fact that 
the red curve in Fig. (|4| is higher than the green curve, and 
this strategy will sort the qubits to be as close the the red line 
as possible. 

As the size of the quantum computer grows, so too will the 
number of entangling operations being carried out in paral- 
lel. Applying the proposed strategy in this Umit of many par- 
allel operations approaches a situation where all vertices are 
paired with other vertices of an identical tilting angle. In this 
limit the expression for the tilting angle after double herald- 
ing (Eqn. ( fTSI l) reduces to the result for when pure, untilted, 
graph are used. Therefore, the average fidelity of GHZ states 
will not deteriorate past the fidelity distribution of building 2- 
qubit tilted graphs. However, the scheme does not mask the 
reduced probability of success. 



E. The reaUgnment procedure 

After a supply of entangled qubits of the required size has 
been produced, those GHZ states that do not meet a criterion 
of acceptable fidelity can be purified by a probabilistic proce- 
dure described in this section. A graphical description of this 
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FIG. 5: The realignment procedure applied to: (a) any qubit in a GHZ state, as all qubits in a GHZ state are cherries; (b) the cherry of the 
inter-node tilted vertex. The rotation required before a computational basis measurement is A/(6'/3), and this is defined by Eqn. l l27t . The 
procedure secedes with probability Ps(/3), Eqn. i28t . Upon failure the tilting in exacerbated, such that the vertices is tilted by angle 
Eqn. (Hill. 



procedure is shown in Fig. (|5]l. One qubit of each low fidelity 
GHZ state is rotated by the unitary matrix M{9): 



M{0) 



-cos{d) sm{d) 
sm{e) coa{e) 



(27) 



For realignment the variable 9 is set equal to 9a, where for 
the rotated qubit cos{9a) is the amplitude of the |0) compo- 
nent. Note that in the ideal limit, as 9a —tt/4, M{9a) be- 
comes the Hadamard gate. Following this rotation the qubit is 
measured in the computational basis. With probability Ps{9a), 
where: 



Ps{9) ^ ^sm\29), 



(28) 



the state |1) is measured, and the remaining qubits are pro- 
jected into a maximally entangled GHZ state. Notice that the 
probability of success is proportional to F^. The expectation 
value E{F^) was proposed in section HTCl as a sensible mea- 
sure of how useful we can expect the product of a given dou- 
ble heralding application. It should now be clear that is a 
good measure, since it tells us how many untilted GHZ states 
we can expect if we attempt to double herald and then attempt 
to realign. 

If the procedure fails then the amount of lopsidedness of 
the GHZ state increases, such that the tilting angle changes to 
—R{9a), where the function is defined such that: 



cos = 



COS^(0) 



1 - isin2(2(/)) 



(29) 



The success probability has an upper bound of 1/2, which 
is approached as cos'^{9a) 1/2. On first inspection this 
contradicts the fact that in the ideal limit an X basis mea- 
surement will deterministically remove one of the qubits to 
give another GHZ state. However, when |0) is measured on 



a state for which cos^{9a) — 0, 1/2, 1, the amount of entan- 
glement lost drops to zero (since, —R{9a) — ~9a, which is 
the same state up to a local Z rotation). If the first attempt at 
realignment fails, a single qubit is lost but the procedure can 
be reattempted on the remaining qubits with a lower success 
probability corresponding to the new tilting angle. 

In addition to purifying tilted GHZ states, the realignment 
procedure is applicable to a wide class of tilted graph states. 
A second example is given in Fig. (HJ)), which begins with the 
tilted graph that was generated as the failed outcome of the 
procedure shown in Fig. The tilted vertex has a single 
neighbour, which itself has no other neighbours, which we 
refer to as a cherry. It is the cherry that is measured out in 
the realignment procedure, and in general, this the is graph 
resource required to attempt to realign a tilted vertex. 

As a closing remark it is worth noting that the rotation 
M{9) and its relationship to R{9) go beyond the realignment 
procedure. Indeed, the M{9) rotation plays a pivotal role in 
both the merge and bridge procedures of later section. Fur- 
thermore, the function R{9) will appear in the state descrip- 
tion of failed procedures. 



III. PROCEDURES FOR JOINING RESOURCE STATES 

So far three things have been demonstrated: (/) during 
phase 1: the degree of tilting in the growth of small GHZ 
states can be limited; (//) at the end of phase 1: tilted GHZ 
states can be realigned into pure graph states (///) when join- 
ing GHZ states, (for phase 2 or 3), there is one available cherry 
with which to attempt realignment. The remainder of this 
paper concerns what happens after (///), that is, what can be 
done with the tilted graph state after the failed realignment of 
Fig. ^p). The naive, and inefficient, answer is that the central 
qubit can be measured out and we can try again. However, in 
this section we propose two procedures that probabilistically 
utilise this tilted vertex. The procedures are coined the merge 



9 




FIG. 6: The starting graph of this figure is the result of a successful 
double-heralding application to two qubits, with Hadamards, from 
GHZ states (with one of the inter-node qubits removed by an X- 
basis measurement). The insets indicate how from this graph we can 
deterministically merge or bridge, A and B. (a) To merge: a X-basis 
measurement is performed on the qubit C; (b) To bridge: a V-basis 
measurement is performed on qubit C. 



and bridge procedure which, when successful, fuse the nodes 
to which the central qubit was connected. Which procedure 
should be used depends on the target graph as both generate 
different kinds of fusion. The effect of a successful proce- 
dure is best seen by comparison with the effect of an X or 
Y measurement on an ideal graph, as shown in Fig. (|6]l. A 
successful bridge procedure will, like a Y-basis measurement, 
connect the nodes by a graph edge. A successful merge pro- 
cedure will, like an X-basis measurement, redundantly encode 
the two nodes in the same way as type-II fusion |30]. 

As a consequence of a failed attempt at either of these pro- 
cedures the resulting state is not a pure graph state. However, 
the failure outcomes can still be described by a further ex- 
tension to our graphical language where failed mergers and 
bridges create partial fusions and weighted graph edges, re- 
spectively. In later sections, we show that even these improper 
graph states are of use, as the partial entanglement established 
by a partial fusion or weighted graph edge can be recycled 
to improve the success probability of subsequent attempts at 
merging or bridging. This extended graphical language will 
be introduced with each procedure. 



A. Merging, and partial fusions 

This section describes how the merge procedures works, 
and how both outcomes are represented as partial fusions. We 
begin with the system corresponding to the failure outcome of 

Fig. m- 



where the first qubit is the tilted central vertex that has x and 
y as neighbours. These are themselves part of two other graph 
segments \ijjx) and In the case of Fig. (|5j)), the graph 

segments have the form: 



1^. 



X 



1 



(31) 



The merging procedure requires only that and \^y) are 
graph states that contain qubits x and y. 

The merger procedure is represented in Fig. and be- 
gins with a rotation of the central qubit by M{±9a), where 
9a parametrises the tilting of the central vertex. This gives the 
state: 



I*) - VpM\l)P.y{±7T/4)\i;x)\i^Y) 



(32) 



where 7 is related to the failure function introduced in the 
realignment section, 7 = ^R{9a), and the operator P^y is 
given by 



Pxyi(t)) = cos((/))l + sin(0)Zj;Zy 



(33) 



For (p = ±7r/4 this reduces to the even (+) or odd (-) parity 
projector. 



1 ± Z.,Zy 



V2 



(34) 



I*) 



icos{9a)\0) + smi9a)\l)Z,Zy) liJxMv), (30) 



Note that we are using the expression parity projector loosely, 
as P^y(±7r/4) is only equal to Pa;j,(±7r/4) up to a constant. 
The reason for the unconventional normalization of Pxy{9) is 
so that when it acts on the graph state \'^y), it generates 
a normalized state, hence the probability of success can be 
read of as the square of the amplitude. 

A measurement in the Z-basis yields a successful merge 
when the |1) state is measured, which happens with prob- 
ability Ps{9). The result is that x and y are projected by 
an even (+) or odd (-) parity projector. Since the choice of 
sign in A/ (±6*) determines the parity subspace of a successful 
outcome, we say that M{±9) targets a particular subspace. 
In the ideal limit where 9 7r/4, the success probability 
Ps{9) 1/2. However, we should expect it to tend to 1 
as the ideal case is deterministic. This paradox is resolved 
when we notice that, as 9 ^ tt / 4, the failure outcome Pxyil) 
tends towards the non-targeted parity projector P^j,(=F7r/4), 
and hence even failure becomes a success. Indeed, the merg- 
ing procedure is completely continuous with an X-basis mea- 
surement, as the prescribed rotation becomes the Hadamard 
under ideal conditions. 

Fig. (|7]i casts this procedure into a graphical language by 
defining a new kind of graph edge, represented by a dashed 
line, that is labelled by an angle 9. We call this dashed line 
a partial fusion, which is to be interpreted as meaning that 
qubits it connects have the operation Pxy{9) applied to them. 
Furthermore, since partial fusions with 9 = ±7r/4 generate 
pure graph states, they are called pure fusions. To see how a 
pure graph state can regained after an impure fusion has oc- 
curred, it is first necessary to see how partial fusions combine. 
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(a) 



Partial fusion notation: 



Apply 

\^ Q to qubits AhB. 



A pure graph state: 

Ha A 



Combing partial fusions: 



008(6*) 



X 




X = COS (6*1 -6*2) 
y = sin(6'i + 6*2) 



H sign(± )=sign(7) 

failure -y \M{±e^\ ^ 
raimre 7^ success 

H» j »' V •!! H» > »-----,-.-« ^ •H 



FIG. 7: The graphical notation for a partial fusion is defined and the procedure for merging tilted vertices is shown, (a) The symbol for a 
partial fusion, a dashed graph edge labelled with an angle. We also show how a maximally entangling partial fusion (Q = ±7r/4) reduces to 
a pure graph state, and how multiple partial fusions combine into a single partial fusion, (b) The procedure for merging tilted vertices when 
there is no pre-existing partial fusion, (c) The procedure for merging tilted vertices when there is a pre-existing partial fusion. The rotation 
required before measuring the tilted vertex is the same in both (b) and (c); M(±S), as in Eqn. \TI\ . Upon failure a partial fusion is generated of 
angle =pi?(0), where K{Q) is the failure function of Eqn. ([29}. Note that, to maximise the probability of success the sign choice in the rotation 
M(±7r/4) must match the sign of the pre-existing partial fusion. 



Consider two partial fusions Pxij{di) and Pxy{d2), which are 
acting on the same qubits as in Fig. (|7^), then: 

Pxy{(l>l)Pxy{(l)2) = COs(0i - (/)2)ll + 8111(01 + 02)^i:^y, 

(35) 

which can be expressed as a new partial fusion ( plus a renor- 
malization constant ): 

Pxy{^l)Pxy{(t>2) - N M{(t>U <t>2)Pxy{(t>) , (36) 

where, 

sin(0) = sin((/)i + (j)2)/NM{(j)iA2), (37) 

and, 

NIM1A2) = 008^(01 -4)2)+ sin^(0i -f 02). (38) 

One significance of these relations is that a pure partial fu- 
sion always overrides a partial fusion, Pxy{±TT / '\)Pxy{(t>) = 
Pxy{±TT/4:), as would be expected. Furthermore, the be- 
haviour of the normalisation constant is crucial since the nor- 
malisation constant determines the probability of success. The 
nature of this dependence will be made explicit in the follow- 
ing disscussion. 

We now consider what can be made of the improper graph 
state which forms the failure outcome of Fig. and the 
procedure that we shall derive is show in Fig. (|7j;). In prepa- 
ration, another attempt at double heralding has to be made on 



two qubits from and When successful, this pro- 

duces another tilted central vertex connected to qubits x and 
y. The state of this system is now: 

I*) = {cosietM + sm{eb)\l)ZxZy) Pxy{-fiMx)\i'Y), 

(39) 

which is the similar to Eqn. dSOb . except for the important ad- 
dition of a partial fusion generated from previous attempts at 
merging, Pxyi'y)- Again, we perform a rotation M{±9b) on 
the central qubit, but this time allow for a choice in sign, the 
purpose of which will become evident soon. This gives the 
state: 

I*) = VPs{0b)\l)Pxy{±7r/4)PxyijlMxMY) (40) 
+ ~ Ps{0bmPxy{T-/2)Pxy{ll)\4'x)\^^Y) , 

which can be simplified using the rules for combing partial 
fusions to give: 

I*) = VpU0^\l)Pxyi±7T/4)\ijx)\^Y) (41) 
+ Vl - Pm{0b,-/imPxy{l3MxMY). 

Here the value for 73 follows from Eqn. JjTl l. and the new 
amplitudes are products of the old amplitudes and the normal- 
ization constant, such that 

Pm{0bni)^Ps{0b)Nlj{±7T/4,Jl) . (42) 
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This simplifies to 

Prni0b,li)=Ps{0b){l±sm{2ji)). (43) 

The second factor of this probability is due to the partial fu- 
sion, and by the correct choice of the sign ± can always be 
made greater than 1 . This is achieved by matching the sign in 
the rotation M{±6b) to the sign of 71. Insight into the phys- 
ical underpinnings of this sign matching can be gained from 
considering that partial fusions can be decomposed into linear 
sums of the odd and even parity projectors. For a partial fusion 
Pxy{4>) we find the following: it is more even parity than odd 
when — 7r/2 < (j) < 0; and more odd parity than even when 
< < 7r/2. As for a rotation M{±9), this will attempt 
to project onto the odd or even parity subspace for + and — 
respectively. Since two graph state qubits initially have equal 
magnitude in the odd/even parity subspaces, a partial fusion 
will increase the magnitude in one particular subspace, and a 
measurement of parity is more likely to work for the dominant 
subspace. 

As a closing remark on merging, notice that if a central 
vertex is not tilted it can be used deterministically to project 
the two qubits into a definite parity state. If the two qubits 
are akeady partially fused then this alters the probability of 
an odd or even parity projection. Some fine tunings of the 
strategy are described in the section ( 1111 Cb . 

B. Bridging, and weighted grapli edges 

When performing a bridge operation on a pure graph state, 
as in Fig. (O, it differs from a merge operation by an addi- 
tional rotation S, where is a diagonal matrix with entries 
(1, z). Again, we shall algebraically describe the bridge proce- 
dure in parallel with a graphical description given in Fig. ([8]), 
and shall begin with a system resulting from the failed out- 
come of Fig. (I5J3). Initially the state is described by Eqn. dSOl l 
and after rotating the central qubit by M{±6a) ■ S the state is: 

l^-) = ^Ps{OamU^y{±^/^)\^x)\^Y) (44) 
+ ^l-ps{eamUxy{lMx)\i^Y) ■ 

Again, 7 is determined by the failure function such that 7 = 
^R{9a), and the new operator is: 

U{(f)) = cos{(f>)i +isin{(j))Z^Zy, (45) 

which is always unitary. 

As with the merging procedure there is a Ps{da) chance of 
successfully obtaining |1) in a measurement, and this gener- 
ates a unitary Uxy{±TT/4:) acting on \iPx)\'4'y)- This unitary 
matrix is identical to a control- Z operation with an additional 
SxSy byproduct. Again, failure establishes a partial amount 
of entanglement that can be recycled in later attempts at bridg- 
ing. 

As with the merging procedure, a failure on the first at- 
tempt is at least partially successful in that U (7) generates 
a control-Z(— 47) gate with a byproduct Zx{—2j)Zy{—2j). 
The notation Z{ip) denotes the diagonal matrix with elements 



(1, exp**^). As before, this partial entanglement will improve 
the probability of success when a second attempt is made. 
Furthermore, as with the merging procedure M^y (0) ■ S and 
Mxy {—0)-S target two distinct but equally acceptable results, 
which generate either a control-Z(7r) or a control-Z(— tt). 

The above procedure is represented in Fig. ((S]?), where the 
partial entanglement of a failed bridge is represented as a 
weighted edge, a solid line labelled with an angle. From the 
definition of Uxy{0), we can again derive a combination rela- 
tion, U xy{di)U xy{02) = Uxy{di + ^2 ) ■ Thcrcforc, unlike the 
partial fusion operator, Uxy (±7r/4) does not override previous 
failures. Consequently, when later attempts are made at bridg- 
ing it is necessary to target the correct amount of weighted 
edge. If there is a preexisting edge of angle 9 then a success- 
ful outcome is generated by 7r/4 — 6* or 7r/4 + 0. As before, 
the sign choice will play an important role in recycling entan- 
glement to boost to the probability of success. 

Taking a system with a weighted edge between x and y and 
then making, via nonideal double heralding, a tilted central 
vertex between x and y, we have the state: 

1^-) = {cosi0b)\O)+smi9i,)\l)ZxZy)Uxy{-fiMx)\4'Y), 

(46) 

A weighted edge of magnitude ±7r/4 — 71 can be targeted by 
rotating the central qubit by M{P±) ■ S, where P± satisfies: 

cos(/3±) = A^B(7i,6'&)cos(6l,,)(±cos(7i)-sin(7i)) . (47) 

Here 

NB{ji,Ob) - (It sin(27i) cos(2(?b))-5 (48) 

plays a analogous role to Nm- After the rotation, the state 
becomes: 

l*> = VPbili,0b)\l)U{±TT/4)\^x)HY) - (49) 
-Ph(7i,^fc)|0)C/(7i + 12)\MHy), 
where Pbiji ,db) is the probability of success 

Pbhi,Ob)=NUj,,9b)ps{0b) . (50) 
The angle 72 of the additional weighted edge equals 

rr^( r. ^^ cos^(6'fc) (±cos(7i) - sin(7i)) 
cos(F(7i,/ii±)) = , , (51) 

^N^\-i^,db)'Ps{()b) 

where -F(7i,/3±) is a generalisation of the failure function 
and it is simple to show that = ^"(0, 0). 
The effect of the previous weighted edge modifies the suc- 
cess probabiUty by a factor Nsi'ji, db)- This factor can al- 
ways be made greater than 1 by the correct choice of the "=f" 
sign, corresponding to a physical choice between in the tar- 
geted unitary [/(±7r/4 — 71). The extent to which Nb devi- 
ates from 1 is dependent on sin(27i), which is intuitive, as it 
vanishes for no pre-existing weighted edge (71 = 0), and it 
is at a maximum for a pure graph edge (71 = 7r/4). The ef- 
fect of 6b on Nb is less intuitive, as it increases as the vertex 
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(a) 



Weighted edge notation: 



Apply U^y{e) 
to qubits X Sz y. 



A pure graph state: 



(b) 



failure 



M{0a).Sr7< 



success 



Combing weighted edges: 



0: 1 



= 6>i + 6>2 



(C) 



failure Ta 





success 



H ±7r/4 H 



FIG. 8: The graphical notation for a weighted edge is defined and the procedure for bridging tilted vertices is shown, (a) The symbol for a 
weighted edge, a graph edge labelled with an angle, is defined. We also show how a maximally entangling partial fusion (6 = ±7r/4) reduces 
to a pure graph state, and how multiple weighted edges combine into a single weighted edge, (b) The procedure for bridging tilted vertices 
when there is no pre-existing weighted. The rotation required before measuring the tilted vertex is M[6a)-S. Upon failure, a weighted edge 
of angle ^R{9a) is created, (c) The procedure for bridging tilted vertices when there is a pre-existing weighted. The required rotation prior to 
measuring the tilted vertex is A/(/3±), where j3± is defined in Eqn. J47| (. Upon failure a weighted edge is generated of angle F(7i, /3±), which 
is defined in Eqn. J51t . Note that, to maximise the probability of success we must choose the correct sign in /3pm. 



approaches a qubit initialization in the |0) or |1) state. Recall 
that a constructively defined graph state generates entangle- 
ment by control-Z gates on initialized qubits, so any qubits in 
the |0) or |1) would not become entangled with the graph. 

Also, we find that the determinism of the ideal regime is 
regained when the central vertex is untilted {0\j — tt/A) even 
if there is a pre-existing weighted edge (71 7^ 0). Although 
the success probability does not become 1, the failure out- 
come becomes the equally desirable alternative to the tar- 
geted operation, that is UxyiTT^/^: — 71) instead of the tar- 
geted UxyiTT^/^: — 71). However, unlike merging in the ideal 
regime, 71 does not alter the probability of these outcomes, 
with each remaining equally likely. 



C. overall strategies 

For both the merge and bridge procedure we have described 
the success probability when using a tilted central vertex of 
angle 6. However, there are alternative methods of generating 
this graph, which would result in a different values of 9. In 
either case, an attempt at double heralding has to be made, 
which will generate a tilted central vertex with a cherry, that 
is tilted by some amount 9a- 

For method (/) the cherry is measured in the Z-basis, so that 
the tilted vertex can be used to merge or bridge with success 



probability: 

P, = N{j,9a)pM, (52) 

where N is the appropriate factor, Nm or Nf,. 

For method (//) an attempt at realignment is made on the 
cherry, which if successful, guarantees success at merging or 
bridging, but if unsuccessful, makes changes the vertex tilting 
to 9a = —R{9a)- Hence, the overall probability of success is: 

Pu ^ I - [I - N{-i,9a)Ps{eam - Psida)], (53) 

which although benefiting from the additional opportunity at 
realignment, may suffer if N{'y,9a)ps{0a) is substantially 
lower than iV(7, 9a)ps{9a)- 

When 7 7^ the best method can be determined by simple 
calculation. On a first attempt at merging or bridging 7 = 
0, and the best method is always the same. If 7 = 0, then 
iV(0, 9) = 1, and since < Psi^a), it follows that P, < 
Pa . Note that on a first attempt Pa has at upper bound of 75%, 
which is approached as 9a 7r/4, whereas on later attempts 
the upper bound can approach 1 . 

IV. IMPROVEMENTS IN RESOURCE COSTS 

The purpose of this paper has been to demonstrate that if 
we know what errors are caused in our entangling operation 
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FIG. 9: The fidelity distribution for the example pair of cavities out- 
lined in appendix |B] The inset is an enlargement of the region just 
below 0.5, with the region of 0.5 - 10"* < F < 0.5 highhghted, 
as this is the window of graph state that are post-selected as being 
within error correction tolerances. 



(due to the additional information provided by detector 'click' 
times) then there are numerous strategies and procedures that 
we can deploy which allow us to construct a graph. This 
scheme is not designed to replace traditional error correction 
(as this will still be necessary), but rather to supplement it. If 
we consider even a modest amount of cavity mismatch, the 
error rate for tilting errors is far beyond the current error cor- 
rection thresholds. Since these errors have to be prevented, the 
only alternative proposal that we are aware of is to post-select 
double heralding applications that meet some fidelity require- 
ment. We give some examples of how our proposal improves 
upon a naive post-selection strategy. 

As a measure of the two approaches we take the prob- 
ability of successfully merging or bridging, on the first at- 
tempt. The pair of cavities we use in our calculations are 
described in appendix |Bl and their photon leakage rates are 
shown in Fig. Taking the fidelity threshold such that 

i — 10"'* < F, the post-selection scheme is constrained 
to the window in Fig. (|9]l. The area under this curve — 
times the rate of photon loss — is the probability of suc- 
cess, which we calculate to be P(post-select) ~ 3.3%. For 
our scheme, we can accept the same graph segments that the 
post-selection strategy accepts, and attempt to merge/bridge 
the lower fidelity graphs, hence P(our-strategy) — P(post- 
select) + P (out- window). On the first attempt at merge/bridge 
Eqn. ( |53] l simplifies to 3P^, and we calculate P(outside- 
window) by integrating 3F^dF over the whole outside region, 
which gives P(outside-window) = 35.7%. Hence, P(our- 
strategy) — 39.0%, and our strategy is an improvement by 
over an order of magnitude. 

Note that, an order of magnitude improvement in gate prob- 
ability equates to far more than an order of magnitude saving 
in the resources (e.g. working memory) required by a quan- 
tum computer Indeed, although the scaling varies between 
different approaches to graph growth, it common to see re- 



source costs that scale with p~ '°s(") where n is a mea- 
sure of the computation size. Hence, an order of magnitude 
improvement in gate success, is equates to an order (0.1)'°^*-"-' 
reduction in resource costs. 



V. CONCLUSION 

In this paper a set of procedures have been proposed that 
allow distributed quantum computing to be achieved despite 
systematic errors that generate a tilting error in the graph ver- 
tices. These errors have been shown to naturally arise when 
the double heralding scheme is used with photon sources that 
possess different inherent leakage rates. Graph construction 
has been considered as consisting in three phases: phase 1 — 
construction of a GHZ resource; phase-2 — efficient merging 
of GHZ resources into larger GHZ states; phase-3 — efficient 
bridging of GHZ states. The proposed scheme contains mod- 
ified protocols for each of these graph construction steps. The 
modified protocols adapt to information from photon detec- 
tion times so as to eliminate all known tilting error We have 
extended the graphical language of graphs states to include 
tilted vertices, weighted edges and partial fusions, so that an 
intuitive grasp of the errors can be gained. Weighted graph 
edges have already found various applications throughout the 
field of quantum computation ll32l l33i [34ll . and the authors 
suspect that tilted vertices and partial fusion may also prove 
to be useful concepts in other problems. 
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APPENDIX A: EVALUATING THE INTEGRAL TO 
DETERMINE 

The expression for E{F'^), the expectation value of the fi- 
delity (minus i) squared, is a complicated function of ti and 
t2 that does not directly integrate to a closed form. However, 
we show that a solution can be found by restricting the ex- 
pression to a region Rj of da and 9^. Furthermore, a solution 
can be also be found for a second region Rj. Rj and Rj 
overlap, and their union covers the entire region of interest, 
< 9a < Tr/2 and < dt < 7r/2. Hence, we can construct 
an analytic solution. First, we express F^ in the following 
form: 



P(p2 



4 



UV 



u + v 

where 8^ are defined in Eqn, 



l + K- 



V 



U + V 



, and: 
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integration. However, first we state the expansion for the sec- 
^ ond region Rj ( tan^iOa) < 2tan2(6'b)): 
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FIG. 10: The analytic 1°* order solution of E{F'^), compared to the 
full exact solution for a specific distribution. The 1^* order solution 
is — up to a factor — independent of the photon leakage rates of the 
cavities. These approximations are represented by solid lines, with 
the colour distinguishing cross-sections of E{F^) along sin^(Sa) ~ 
sm^{6b) and sm^{Oa) = 1 — am^{9i,). For the first cross-section, 
in red, the solution is exact even up to 0*'' order For the second 
cross-section, in green, we compare the approximation with the a 
full solution, which is represented by a dashed line. For the full 
solution we use the same cavity parameters as in Fig. (|2};,|4l, which 
is specified in appendix IB] 



and the functions U and V are: 



U = PA{tl)PB{t2), 
V = PB{tl)PA{t2). 



(A3) 



Since U and V are never negative, V/ {V + U) is always less 
than one. Therefore, provided /-iT < 1 we will be able to per- 
form a binomial expansion. The restriction on K gives us a 
region Rj, ta.v?{6b) < 2 tan^(0a). The expansion is: 



n—1 ^ ^ 



E{F')R,=J2{-irQ2i-^-l] Jn, (A6) 

n—1 \ 1 / 



where, 



dtidt2 



(A7) 



Since the difference between U and V is only in the labelHng 
of ti and t2, we know J^=I^. Furthermore: 
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Interestingly, for odd n. In are fixed in terms of Ik with even 
k. However, the even terms of /„ are not constrained by this 
recurrence relation. For example, Ii = i/g, which means that 
(up to a constant factor Iq) we can expand up to first order 
without performing any numerical integration. Hence, when 



sm 



9,)-sin^(0b)«O, 



/o, (A9) 



where the change between E{F'^)rj and E{F'^)r,, is 
achieved by using 8l as the largest of the pair 9i & 02, 
and 05 as the smallest. Two cross-sections of E{F^)aprx are 
shown in Fig. ( fTOb . and compared with an example full solu- 
tion. Both cross sections agree well with the exact solution, 
even far away from sin^(6'a) — sin^{6i,) — 0. The importance 
of this result is that it tells us about the shape of E{F'^) inde- 
pendently of the shape of PA{t) and Psit). In the main body 
of this paper we argue from a numerical solution of E{F'^) 
— shown in Fig. (HJi — that sin'^ (6'a) = sin^ {6b) is prefer- 
able to sin^ 9a = I — sm^{9b)- The above result shows us 
that, at least near 9a ~ 7r/4, this behaviour is not an arte- 
fact of the particular solution but a general property valid for 
any PA{t) and Psit). In fact, the only requirement is that Pi 
are non-negative, so we can also accomodate photon loss and 
distributions that violate our Hamiltonian, Eqn. Q. 



APPENDIX B: EXAMPLE CAVITY SOLUTIONS 



where. 

An expression for E{F'^)fjj can be found by calculating 
/„. Although these may have to be performed numerically, 
this is preferable to numerically calculating E{F'^) at every 
point. Also, we shall show that there is a recurrence relation 
between different /„, so some terms can be found without any 



Here we outline solutions to the Jaynes-Cummings model 
and give an example solution that is used throughout the pa- 
per for expository purposes. Solutions to Eqn. (|3]l, involve 
two parameters, gx and Hx- In the strong coupling regime 
{Kx ^ '^gx), the system tends to experience Rabi oscillations 
before the photon escapes from the system. These oscillations 
will appear in the probability distribution for photon detection. 
On the other hand, in the weak coupling regime (k^; 3> 4g^), 
the system experiences the quantum Zeno effect due to con- 
tinual measurement of the cavity, and the relaxation time for 
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the atom-cavity system increases with as t ^ n^j g\- Ar- 
guably, this effect is an artefact of the model, as when the 
coupling decreases, it only decreases with respect to the cav- 
ity mode. Therefore, in the weak coupling regime the system 
will lose photons faster than predicted, just into modes that 
are not being monitored with detectors. However, this is not a 
major concern for our paper, as the primary purpose of invok- 
ing the Jaynes-Cummings Hamiltonian is to provide a means 
of modelling the collapse process. The actual photon leakage 
rates, Px{t), will be determined experimentally. 

Concerning the example solution: Rabi oscillations are un- 
desirable for a system being used in double heralding; and the 
model is unsuitable for very weak couphng. Between these 
regimes there exists a class of critically damped solutions — 
for which k.^ = ^Qx — which have a convenient form that we 



take advantage of: 

Cx{t,gx) = Px{t,gx,4gx) = Ag^t^ eicp-^^^' Q{t), (Bl) 
where Q{t) is the Heavyside function 

The example pair of cavities used throughout this paper, are 
parametrized by = 10 and gA = 12.5. Note, that all con- 
clusions of this paper are independent of the particular form 
of the photon leakage rate, and that this example is solely used 
to guide the reader through the paper. 
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